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ABSTRACT 


The development and structure of Dean vortices in a high aspect ratio curved 
channel, were measured and studied for Dean numbers from 35 to fully turbulent con¬ 
ditions. Mean vorticity components were determined from measurements of the three 
components of mean velocity, which were measured using a miniature five-hole pressure 
probe. Streamwise, radial and spanwise vorticity components measured 120° from the 
start of curvature, show evidence of Dean vortices in the time-averaged flow field. Peak 
streamwise and peak spanwise vorticity, as well as streamwise circulation, generally in¬ 
crease as the Dean number increases from 60 to 240. Peak radial vorticity increases as 
Dean numbers vary from 60 to about 160. Spectra obtained from hot-wire probes show 
evidence of different frequency events, including twisting and undulation. Twisting is 
apparent for Dean numbers of 135 to 160 while undulation is observed at lower Dean 
numbers. Flow visualization results illustrate a variety of time dependent events in the 
channel. 
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I. INTRODUCTION 


Fully developed laminar channel flow does not maintain a parabolic velocity profile. 
Centrifugal instabilities in the flow result in counter-rotatw - vortex pairs, which initially 
form near the concave surface. The imposition of the centrifugal instabilities begin be¬ 
cause centrifugal forces are greater on the fluid particles at the center of the channel than 
near the walls. Consequently, faster moving particles near the center are forced to move 
toward the concave surface, displacing the fluid particles located there. These displaced 
fluid particles move spanwise and then away from the concave surface. The path they 
follow is then like an arc, as shown in Figures 1 and 2. The resulting counter-rotating 
vortex pairs are known as Dean vortex pairs. 

Even though flows near compressor blades and turbine blades are not channel flows, 
centrifugal instabilities result in the formation of counter-rotating vortex pairs near 
concave surfaces in all three situations. Similar effects are also present in combustion 
chambers, rocket nozzles as well as internal cooling passages of turbine blades. Pre¬ 
dicting curved channel flow is an attractive first step in accounting for curvature in the 
practical environments because channel flow is a much simpler undertaking. Develop¬ 
ment of means to predict flows near concave surfaces in the practical environments will 
eventually lead to improved methods of protecting surfaces from exposure to hot gases. 
This will eventually allow higher inlet temperatures which result in greater efficiencies, 
improved power-to-weight ratios, and better engineering designs. 

A. BACKGROUND 

The first analytical study of curved channel {d4.r,) flow was performed by W.R. 
Dean [Ref. 1] in 1928. Dean found that the flow becomes unstable due to the effect of 
concave curvature when small disturbances are present. He also found that the critical 


I 





minimum flow velocity at which the small disturbances start to grow occurs when 
Ud/\\'djr, is above 36 for a channel with radius ratio near 1.0. This parameter is now' 
know'n as the Dean number, De, where U is the mean velocity in the channel, v is the 
kinematic viscosity, d is the spacing the between inner and outer walls of the channel, 
and r , is the radius of curvature of the inner wall. Reid [Ref. 2] provides additional in¬ 
formation on the neutral stability curve. Using a low aspect ratio curved channel, 
Hawthorne [Ref. 3] describes streamw'ise velocity distributions which result in curved 
channel flow'. Cheng et al. [Ref. 4] show how' flow’ visualization patterns in 
radial spanwise planes change w'ith aspect ratios, where aspect ratio is varied from 1 to 
12. Works by McKee [Ref 5] and Flentie [Ref. 6] resulted in the first published work 
on high aspect ratio curved channel flow by Kelleher et al. [Ref 7], They show' how' 
streamwise velocity and Dean vortex pair spacing varies with Dean number. Finlay 
et al. [Refs. 8, 9 ] present a comprehensive analytical and numerical investigation of 
Dean vortices in high aspect ratio curved channel flow'. Their work predicts two kinds 
of wavy Dean vortex flows, referred to as undulating Dean vortex flow' and twisting 
Dean vortex flow [Ref. 10]. Experiments by Ligrani and Niver [Ref 11], Finlay, Ligrani 
and Bland [Ref 10], and Alfredsson and Persson [Ref. 12] confirm the presence of these 
secondary instabilities. Ligrani et al. [Refs. 11,13], Niver [Ref 14], Longest [Ref. 15], 
Ligrani et al. [Ref. 16], and Baun [Ref 17] present additional results from the same 
curved channel used in the present study which has an aspect ratio of 40 to 1. A recent 
w'ork by Finlay, Ligrani and Bland [Ref 10] compares the numerical and experimental 
results when wavy vortex flow is present. Ligrani and Longest [Ref. 13] describe vortex 
pair appearance and disappearance events in curved channel flow. 

B. OBJECTIVES 

The primary objective of the present study is to measure flow characteristics in a 
curved channel with 40 to 1 aspect ratio at Dean numbers ranging from 35 to 425. A 
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miniature five-hole probe is used to obtain surveys of a number of time-averaged quan¬ 
tities, including total pressure, streamwise velocity, secondary flow vectors, streamwise 
vorticity, radial vorticity, and spanwise vorticity. Power spectra of the voltage fluctu¬ 
ations from a hot-wire probe sensing velocity fluctuations are used to show the fre¬ 
quencies of ■. iflerent events in the flow field. Unsteady Dean vortex behavior is 
illustrated bv flow visualization results obtained over radial/spanwise planes at three 
streamwise locations. 

C. ORGANIZATION 

Chapter II describes curved channel. Chapter III describes the equipment and 
procedures used for measurement of pressure and velocity surveys. Chapter IV describes 
the procedures and equipment used to obtain power spectra. Chapter V gives equipment 
and procedures used to acquire flow visualization results. Experimental results are pre¬ 
sented in Chapter VI. Summary and conclusions are presented in Chapter VII. Ap¬ 
pendix A contains all of the figures. Appendix B presents the results of an uncertainty 
analysis from Baun [Ref. 17|. A software directory of the programs used in this study 
are presented in Appendix C. 
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II. CURVED CHANNEL 


A. CHANNEL DESCRIPTION 

The transparent curved channel is located in the laboratories of the Mechanical 
Engineering Department at the Naval Postgraduate School. Detailed explanations of 
the design and construction of the facility are found in Sieband [Ref. 18], Niver [Ref. 
14], Ligrani and Niver [Ref. 11], and Baun [Ref.17], A schematic illustration [Ref. 14] 
and photograph [Ref. 18] of the rectangular curved channel are shown in Figures 3 and 
4, respectively. The test facility equipment between plenum #1 and the blower exhaust 
is correctly represented in Figure 3. The photograph, Figure 4, shows a plenum/blower 
configuration that is no longer used. For the sake of continuity, a brief discussion of the 
channel is provided below. 

The channel is an open circuit suction facility designed for low speed transition 
studies. Its primary function is in the study of the structural characteristics of the flow 
e.g., pressure, velocity, vorticitv distributions, and flow visualization results. It is trans¬ 
parent because it is constructed of Lexan, a polycarbonate. The lip at the channel inlet 
is constructed of quarter circumference sections of 15.2 cm outside diameter pipe. To 
reduce spatial non-uniformities in the flow, an aluminum honeycomb and three fine 
mesh screens follow. This section is followed by a 20:1 contraction nozzle, whose con¬ 
tour is given by a fifth order polynomial. A photograph of the inlet section and nozzle 
is seen in Figure 5. The flow is accelerated through the nozzle into a 2.44 m long 
straight duct. This entry length provides fully developed flow at the inlet side of the 
curved section for Dean numbers less than 640 [Ref. 18]. The curved test section has 
an inner to outer radius ratio of 0.979 and turns the flow 180 degrees. The flow then 
enters another 2.44 m straight section before reaching another honeycomb, screens, and 
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a dilfuser. The diffuser is followed by the outlet and blower plenums A gate valve and 
orifice plate are used to control and regulate the flow in the channel, and are mounted 
between these two plenums on a 5.25 cm ID (2 inch schedule 40) pipe. A low pressure 
is maintained in the blower plenum by a 250W (one-third horsepower) motor. Rubber 
mounts and flexible couplings between the plenum and the blower as well as beneath the 
blower, isolate the channel from mechanical vibration. The side walls, longitudinal and 
cross beam supports maintain a cross section of 1.27 + 0.015 cm by 50.8 ± 0.05 cm. 
[Ref. 16] 

Flow' properties in a radial'spanwise plane are measured after a miniature five-hole 
pressure probe is inserted through a slot cut in the curved section. The slot is 0.32 cm 
wide and 7.62 cm long and is located in the convex w'all of the channel. The slot is 
aligned in the spanwise direction 5.08 cm from the centerline at a sheamwise location 
of 6 = 120°. A support block made of Lexan maintains the convex wall structural in¬ 
tegrity and minimizes the deflection of the curved wall in the vicinity of the slot. The 
slot is lined with foam to eliminate air leakage when the pressure probe is inserted. A 
schematic illustration of the Lexan support block and foam sealing system [Ref. 17] is 
shown in Figure 6. 

B. DEAN NUMBER DETERMINATION 

The mass llow r rate and Dean number are determined from measurements of the 
pressure drop across the AS ME orifice plate. The pressure drop is measured with a 
Validyne Model PS309 digital pressure manometer. The signal is collected and processed 
using a Hewlett-Packard (HP) 3498A extender, HP 3497A data acquisition/control unit 
and a HP 300 Model 9000 computer. Additional details on the data acquisition system 
and computer are given in the following chapter. 

The pressure drop is used to make a first guess at the Dean number by interpolating 
pressure drop versus Dean number data from Niver [Ref. 14]. Density is calculated using 
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the ideal gas law and ambient conditions. The expansion coefficient, Y, is calculated 
from Holman and Gajda [Ref. 19]. An initial mass flow rate is calculated from the initial 
estimate of the Dean number and pipe Reynolds number, Re p . The latter is calculated 
using: 


rnd 

Re p = -f 

p pvA, 


{ 1 } 


The flow coefficient, K, is found by interpolation from the ASME tables [Ref. 20 ]. The 
mass flow rate is then updated using: 

m = KA 0T Yj2pAP 0T {2} 


Iterations of this procedure are continued until the mass flow rate converges to within 
one percent of the previous value. A final estimate of the Dean number is then calcu¬ 
lated using: 


De = 



{3} 


To set the curved channel up for a run at a particular Dean number, this procedure is 
implemented using the computer program DEAN 15 described in the Appendix C. [Ref. 
17] 
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III. PROCEDURES FOR MEASUREMENT OF PRESSURE AND 


VELOCITY SURVEYS 

A detailed description of the fabrication, use, and calibration of the five-hole pres¬ 
sure probe is given by Ligrani et ah [Refs. 13,16, 21 }. For completeness in thought and 
context, portions of these articles are presented below. 

A. PROBE DESCRIPTIONS 

1. KIEL PROBE 

A miniature Kiel probe is used to measure time-averaged total pressure distrib¬ 
utions, in a radial, spanwise plane, 120° from the start of curvature for Dean numbers 
from 50 to 450. This data is taken to qualify flow behavior of the channel after a thor¬ 
ough cleaning and reworking of the inlet apparatus. Secondly, these data are used to 
validate the performance of the five-hole pressure probe. The diameters of the Kiel 
probe shield and pressure port are 1.68 mm and 0.4 mm, respectively. 

2. FIVE-HOLE PRESSURE PROBE 

Two schematic illustrations and a photograph of the miniature five-hole pres¬ 
sure probe [Ref. 17] are shown in Figures 7, 8 and 9, respectively. In each, five tubes 
are evident such that the central tube is surrounded by four other tubes. The probe tip 
diameter is 1.38 mm, a value different from Baun [Ref. 17] (1.22 mm) because the probe 
was reconstructed by fitting a new tip to the original mounting cylinder. The end of each 
of the surrounding tubes is beveled at a 45° angle with respect to the centerline of the 
center tube. Each tube is constructed from stainless steel hypodermic tubing with inner 
and outer diameters of 0.203 mm and 0.406 mm. The five tubes were first joined by spot 
welds, bent in a 90° turn with a radius of about 1.5 mm and then silver soldered along 
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their length. Finally, a surface grinder was used to bevel the four outer tubes. This 
approach is slightly different from the one described by Treaster and Houtz [Ref. 22]. 

The probe tip extends 5.08 mm (four tip diameters) beyond the bend, as shown 
in Figure 7. This minimizes the effects of stem flow blockage at the measuring location. 
Similarly, the tip is 22.9 mm (19 tip diameters) from a 9.53 mm diameter mounting cyl¬ 
inder. This mounting cylinder was sized for easy mounting and handling. 

The miniature five-hole pressure probe is designed to minimize the adverse ef¬ 
fects of flow blockage and spatial resolution. When the probe is in place for measure¬ 
ments, the mounting cylinder is not exposed to the flow. A five hole rather than a four 
hole design was chosen because it is easier to calibrate, and because spatial resolution 
corrections are more readily applied (Ref. 21}. The tip is conical rather than prismatic 
to minimize flow blockage effects when the probe is near the wall. 

B. EXPERIMENTAL DETAILS 

Imperial Eastman 22-P-l 8 polyflo tubing is used to connect the pressure probes to 
the pressure tramducers. To avoid pressure oscillations due to gas compressibility and 
flow passage constriction [Refs. 16,17], the maximum length of tubing used is 18 cm. 

For pressure measurements using the Kiel probe, a Celesco Model LCVR variable 
reluctance differential pressure transducer is used with a full-scale range of 2 cm of wa¬ 
ter. Each port of the five-hole probe is connected to a Validyne Model DP103-06 vari¬ 
able reluctance differential pressure transducer with a full-scale range of 0.25 cm of 
water. At Dean numbers greater than 370, pressure ranges exceed those of the Validyne 
transducers and so, Celesco transducers are used instead. Each transducer signal is sent 
to a Celesco Model CD01D carrier demodulator, which converts outputs from the 
transducers to DC voltage signals which are proportional to pressure. 

Using the computer program PROCAL, a calibration procedure similar to that de¬ 
scribed by Baun [Ref. 17] is used, except that a Validyne Model PS309 digital pressure 
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manometer is used instead cf a standard fluid manometer. As expected, the calibration 
coefficients (cm H 2 Ojvoli) o f the Celesco transducers are ten times greater than those of 
the Validyne transducers. 

The voltage signals are r ent through an IIP 3498A extender to an HP 3497A data 
acquisition control unit. These units then send data signals to the HP 300 Model 9000 
computer to be processed and stored. In a procedure similar to that used by Baun [Ref. 
17], 50 voltage signals are collected sequentially from each carrier demodulator. The 
mean and standard deviation of voltages from each channel are calculated. Sampling 
continues until 50 continuous samples from a particular channel are within 2.5 standard 
deviations of the previously calculated mean value. 

C. PROBE CALIBRATION PROCEDURES AND RESULTS 

The miniature five-hole probe is calibrated using procedures closely resembling those 
used by Baun [Ref. 17). This is accomplished in an open circuit, subsonic wind tunnel 
illustrated in Figure 11. Manufactured by ALROLAB, the tunnel is located in the same 
Mechanical Engineering Laboratory as the curved channel. The calibration sled, shown 
in Figure 10, is mounted on the top of the walls of the test section. Before calibration 
begins, the r robe is positioned in the flow using the sled so that yaw angle is zero. The 
computer program. ORIENT is used to ai' J * ins procedure. The probe is considered to 
be at zero pitch angle when positioned vertically in the calibration sled. Calibration is 
undertaken using twenty-five angle combinations, i.e., five yaw angles at each of five 
pitch angles. For each angle combination, yaw, pitch, total, and total minus static 
pressure coefficients are calculated using equations given by [Refs. 16,17]: 


C py = (P 2 -P 3 )l(P ] -P) 

H) 

C pp = (P i -P,)t(P x -P) 

{5} 
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C pl = (P,-P)!(P ] -P) 


{ 6 } 


C p!s = {P,-P s )HP\-P) {7} 

and 

C p ^{P x -P s )i{P,-P) {S} 

Probe calibration results which correspond to a mean velocity of 3.9 m/s are pre¬ 
sented in Figures 13 through 17. Figure 13 shows increasing values of C py with in¬ 
creasing yaw angle, with C py = 0 for zero yaw and pitch angles. Figure 14 shows C pp 
increasing with increasing pitch angle, but shows C py = — 0.25 at zero pitch angle for all 
values of yaw angle. In Figure 15, C p ,, and in Figure 16, C p „ , both show a dependence 
on yaw angle at a given pitch angle. Figure 17 shows C p , increasing slightly with pitch 
angle. 

D. CURVED CHANNEL MEASUREMENT PROCEDURES 
1. TRAVERSING MECHANISM 

The probe is held in place within the channel using a mount which allows initial 
yaw angles and pitch angles to be adjusted within ±0.1° and + 0.5°, respectively. This 
mount is attached to a travelling block whose position is controlled by an automated 
two-dimensionaj traversing mechanism [Ref 17j. A photograph [Ref. 17) of this equip¬ 
ment is seen in Figure 12. The block is moved using two leadsciews with a pitch of 1.27 
mm (0.05 inch). A Superior Electric Type M092-FD-310 stepping,'synchronous motor 
rotates each leadscrew'. These motors are controlled by a Superior Electric MITAS Type 
PMS085-C2AR controller and a Superior Electric MITAS Type PMS085-D050 drive. 
These equipment direct the probe movement as directed by the HP 300 Modei 9000 
computer. 
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2. MEASUREMENT PROCEDURE 


The probe is initially placed in the mount with the travelling block as close to 
the channel curved surface as possible. The probe is manually oriented so that it barely 
touches the concave wall. The probe mount is tightened and the probe is moved pre¬ 
cisely to a radial location of y/d = 0.5 by the M1TAS controller. The computer program 
ALIGN is then used to align the probe to zero yaw and pitch angles. When the probe 
is pioperly aligned with zero yaw angle and zero pitch angle, the pressure readings on 
opposing ports (i.e., 2 and 3, 4 and 5) are equal. Fine adjustments to probe orientation 
are accomplished by a specialized mount (seen in Figure 12), designed by Baun [Ref. 

17]. 

Data collection begins with the probe 0.127 cm from the concave wall at the 
first spanwise position. The probe moves toward the convex wall as measurements are 
made. To move to the next spanwise position, the probe is moved 0.254 cm in the 
spanwise direction, then 0.127 cm in the reverse spanwise direction and finally toward 
the concave wall until it is 0.127 cm away. This positioning procedure is used to reduce 
the effects of leadscrew backlash and ensure repeatable data regardless of the direction 
of spanwise probe movement.[Ref. 17] 

At each probe position in the channel, 100 pressures are sampled per pressure 
port, at a frequency of about five Hz. After the probe has moved to the next location, 
measurements are delayed ten seconds to avoid sampling when unsteadiness caused by 
probe movement is present in the flow. A full survey of pressure measurements using 
the five-hole probe requires about eight hours. Because the Kiel probe involves meas¬ 
urements of only one port, the elapsed time is much less. The entire traversing operation 
is computer controlled and performed late at night to minimize channel inlet disturb¬ 
ances caused by any sort of motion in the building. 
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3. DATA REDUCTION 


Procedures similar to described used by Ligrani ei al. [Refs. 16,21] are used for 
processing of data. The first step involves correcting pressures to account for spatial 
resolution effects using procedures described in [Ref. 21]. C n and C pp are then determined 
using the same procedures described in [Ref. 16], except for the preliminary’ pitch calcu¬ 
lation. For this, the following equation is employed: 

p' = 8.458(C„ + 0.2568) - Poff {9} 


p off is zero for the five-hole probe arrangement used in this present study. After the yaw 
and pitch angles are known, a double interpolation procedure is employed to compute 
C,„ and C pl . Total velocity and three velocity components are then calculated. The final 
phase involves adjustments made to account for pressure and velocity gradients meas¬ 
ured in the three-dimensional flow field (Ref. 21]. 

Using finite difference equations, distributions of the three components of mean 
velocity determine the three components of mean vorticity. Secondary flow vector 
magnitudes are used to calculate streamwise, radial and spanwise vorticity, using the 
following equations [Ref. 10]: 


(Da = 


8U r dU 2 


<d z = r 


cz 

8r 


dU e 


oz 

V e 

dU e 


dr 


{ 10 } 


{ 11 } 


{ 12 } 
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U e , U r and U 2 represent the streamwise, radial and spanwise components of velocity. 
The streamwise derivatives are neglected for the determination of w, and cj,. The first 
term in equation {12}, UJr, is included in the determination of spanwise vorticity but 
was neglected in the study by Baun [Ref. 17]. 
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IV. POWER SPECTRA PROCEDURES AND EQUIPMENT 


Hot-wire anemometers are widely used for measurement of instantaneous velocity 
components in turbulent gas flows. Their operation relies on the fact that the electrical 
resistance of a metal conductor is a function of its temperature [Ref. 23]. The sensing 
element of a hot-wire probe is a very small diameter wire, usually made of platinum 
and/or tungsten. As such, the wire has high spatial resolution and provides little inter¬ 
ference to flow. Hot-wire sensors are heated by electric current, where the amount of 
heat loss to the gas flow' is directly related to the flow velocity. This relationship is de¬ 
termined by calibration [Ref. 24], 

In this study, the hot-wire is used in the constant-temperature (or constant- 
resistance) mode of operation. As such, the hot-wire sensor forms one leg of a 
Wheatstone bridge. During operation, this bridge circuit is balanced with the hot-wire 
exposed to a zero-velocity flow field, and the resistance required to balance the bridge 
is the cold resistance. When the hot-wire is operated at elevated temperatures, its re¬ 
sistance is the hot resistance, determined by multiplying the cold resistance and the 
overheat ratio. Standard values of the overheat ratio range from 1.2 to 1.8. When hot¬ 
wire sensors are exposed to fluctuating velocities, a variable current feedback amplifier 
maintains the balance of the bridge, as well as the sensor resistance, by passing different 
levels of current to heat or cool the wire so that the average temperature along its length 
is maintained constant. The voltage required to keep the sensor resistance constant is 
then used as a measure of flow velocity fluctuations in the vicinity of the sensor. 

The widespread use of hot-wire anemometry systems is due to several advantages. 
First, the small dimensions of the sensor provide very little flow blockage and alteration. 
Secondly, electronic system response is almost instantaneous and rapid fluctuations in 
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velocity components can be recorded with minimal electronic distortion. Third, system 
sensitivity is such that detection of fluctuations of only a few percent of the mean ve¬ 
locity can be recorded. Finally, electrical output of a hot-wire anemometer can be 
r caJi;y piocessed by both digital and analog systems.[Ref. 23] 

A. HOT-WIRE PROBE 

A DAXTEC Type 55P01, conventional hot-wire probe is used for the present study. 
The sensing wire is fabricated of platinum plated tungsten wire, five microns in diameter 
with a sensitive wire length of 1.25 mm. The wire is oriented so as to be sensitive to both 
the streamwise and spanwise velocity components. The probe is held in place by a 
DAXTEC probe mount which is attached to a plug in the channel side wall. The lo¬ 
cation of the probe in the channel is d= 112°, y/d=0.S with a spanwise location of 9.84 
cm (3.875 in) or 10.16 cm (4.0 in) from the channel side wall. A photograph of the 
mounted hot-wire sensor is given in Figure 77. 

B. SIGNAL MEASUREMENT 

During measurements, the hot-wire probe is connected to a DAXTEC 55M01 
Standard Bridge using a five meter coaxial cable. The system is operated in the 
constant-temperature mode with an overheat ratio of 1.8, bridge gain of 3.0, and an HF 
FILTER setting of 2.0. The output signal from the DAXTEC Bridge is passed through 
a DAXTEC Model 56X20 Signal Conditioner, which filters and amplifies the voltage 
signal before it is sent to the signal analyzer and a B/K Precision Model 1476A 10 MHz 
oscilloscope. The high-pass filter of the signal conditioner is set at 0.1 Hz to remove the 
DC signal. To prevent aliasing, the low-pass filter is set at 0.1 or 0.3 kHz, depending 
on the frequency range being investigated. The gain setting is 100. 
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C. SIGNAL PROCESSING AND ANALYSIS 

An HP 3562A Signal Analyzer is used to obtain spectra of signals from the hot-wire 
anemometer. The HP 3562A is a dynamic, dual-channel fast Fourier transformed (FFT) 
based analyzer with a 26.5 ^Hz-to-100 kHz frequency range and a 150 dB measurement 
range. The nominal maximum sampling rate for all channels is 256 kHz. A schematic 
diagram of the measurement process [Ref. 25] is given in Figure 78, and is also described 
below. 

The analog to digital conversion (ADC) portion of the analyzer works in two steps. 
First, it takes a picture of the signal at one point in time. Second, this analog sample is 
convened to 14-bit digital words in a time record which is 2048 samples long. The 
measurement stage then begins after the analog signal is converted to digital data and 
properly scaled. In this part of the measurement process, data are filtered to an appro¬ 
priate frequency span. Windowing is then performed using software and two specialized 
processors. The fast Fourier transform is subsequently implemented to convert the time 
vary ing signal to a frequency varying one. Next, the floating point processor (FPP) then 
provides the computational power necessary for averaging and calculating measurement 
displays. Signal averaging concludes the measurement portion and the display portion 
begins. 

The vector display used in the HP 3562A employs 1601 display points, so data must 
be compressed or expanded to fit on the display, as appropriate. With the display, the 
input data (view input), filtered input data (filtered input), or output data (measurement) 
may be presented. Input data are ones extracted at the output of the HP 3562A 
analog-to-digital converter, which have not been through the digital filter or measure¬ 
ment process. Filtered input data displays show input signals after they have been dig¬ 
itally filtered. Displays of output data show results after signals are filtered and 
processed according to settings selected by the user. [Ref. 26] 
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One channel is employed in the present study. The analyzer is used in the power 
spectrum mode of display with the trigger set to free run. A Hanning data window 
(Hann window) is utilized to attenuate the input signal at both ends of the record, forc¬ 
ing the signal to appear periodic. Time Averaging Off is selected so that the analyzer 
makes continual measurements with one average, erasing the result of each previous 
measurement. Time traces of the input signal are taken from the filtered input displays 
and are plotted in terms of volts versus seconds. Power spectra, taken from the meas¬ 
urement displays, are plotted with axes of (Voits 2 l Hz) versus Hz. Typical settings of the 
HP 3562A Signal Analyzer used in this study are presented in Figure 79. A photograph 
showing the electronic components used to obtain the power spectra is given in 
Figure 80. 
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V. FLOW VISUALIZATION APPARATUS AND PROCEDURE 


The apparatus used for flow visualization are now described, including the smoke 
generator, camera, lighting system, video recorder, video printer, and video monitor. 

A. SMOKE GENERATOR 

The smoke generator, which was first used by Morrison [Ref. 27], is shown sche¬ 
matically in Figure 96 [Ref. 28]. It consists of five main components: combustion 
chamber, reaction vessel, condensers, moisture separators, and smoke rake. Combustion 
takes place inside a 40.64 cm long, 6.27 cm inside diameter (2.5 inch schedule 40) steel 
pipe. Inside this chamber, smoke is produced from a mixture of hickory and mesquite 
wood chips, mixed in approximately equal amounts. The wood is initially ignited by 
passing current through a coil of several strands of 0.03175 cm (0.0125 in) diameter 
nickel-chromium wire, which is located beneath the wood chips, as compressed air at a 
pressure of 2-8 psig is passed through the pipe and between the wood chips. Different 
flow rates of smoke are produced by regulating the pressure level, and hence, the air 
supply to the combustion pipe. At the exit of the pipe, the smoke flows into a 1500 ml 
Pyrex reaction vessel where some unwanted combustion particles are removed. Next, it 
flows through two 44.5 cm Pyrex coiled Graham condensers where the smoke is cooled 
to temperatures near ambient. Two moisture separators, arranged in series, follow the 
condenser and allow additional moisture to be removed from the smoke. Each moisture 
separator consists of a cylindrical cavity 13.34 cm long and 11.43 cm in diameter, with 
loosely packed aluminum wool to increase the mass transfer surface area. These sepa¬ 
rators are connected to 1.22 cm inside diameter rubber tubing which carries the smoke 
to a manifold and rake located at the inlet of the curved channel. The rake contains nine 


18 





0.34 cm inside diameter tubes connected to the manifold which has an inside diameter 
of 2.24 cm. 

B. CAMERA AND LIGHTING 

A Smith-Victor 710-SG photographic floodlight is used to illuminate the smoke as 
flow visualization experiments are conducted. It is positioned so that light is directed 
radially outward toward the convex side of the curved channel. To prevent the Lexan 
of the curved section from deforming due to heat from the floodlight (it is rated at 
600W), 2.54 cm thick, foam insulating material is placed along the outside of the convex 
wall. 

A Dage-MTI CCD-72 Series Video Camera with an AMPEX Vidicon, 1:1.4 f = 25 
mm lens is used to photograph the flow patterns formed by the smoke. The camera head 
is connected to a control unit by means of a cable assembly. The camera is solid state 
and utilizes a high resolution CCD (Close Coupled Device) imager and high performance 
analog processor to provide an excellent combination of precision and flexibility. The 
camera head is quite small measuring 3.05 cm by 4.32 cm by 12.2 cm long and weighs 
180 grams. The Camera Control Unit (CCU) has an IR cut filter and optical low-pass 
filter to minimize aliasing, and provides controls for gain and black level adjustment. 
An illustration of these controls on the front panel of the CCU are evident in 
Figure 97. Also apparent are the means to adjust bandwidth control, enhancement 
control, polarity reversal, gray scale stretch, and gamma settings. A test signal is also 
available to allow' adjustment of the contrast of the image on a monitor [Ref. 29]. 
During this study, the gain is set at 4.62, the black level is set at 8.88, and the polarity 
switch is set to the positive position. Because the gain control is used in manual, the 
gray scale stretch feature is not used. For optimum sharpness with minimum noise, the 
bandwidth and edge enhancement settings are 2/3 of the maximum setting, and the 
gamma correction is set at 0.70. 
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A schematic showing the camera and lighting setup is presented in Figure 95 [Ref. 
17]. With this arrangement, photographs of smoke patterns in radial/spanwise planes 
are taken at three locations on the curved section: 0 = 40°, 80° , and 120°. 

C. VIDEO RECORDER, PRINTER AND MONITOR 

In order to obtain time histories of flow patterns, the CCU is connected to a 
Panasonic AG-1960 S-VHS video cassette recorder, which operates in the Super-VHS 
format, using a TDK XP-ST120 videotape. 7 cm by 10 cm hard copies of the individual 
images from the video are produced using a Mitsubishi P65U Video Copy Processor 
(VCP). The VCP allows the user to control contrast, print orientation, and scan type, 
as well as black and white levels, picture sharpness, gain and filter levels. In this study, 
the contrast and print orientation switches are set to normal. Scan selection is set to the 
field mode to best capture the rapidly varying flow phenomena. The black and white 
levels are varied to compensate for lighting differences due to smoke density. The 
sharpness, gain and filter selection are set to off Print size and position may also be 
controlled using an additional menu. Here, the standard (default) settings for print size 
and position are utilized. 

The video image is monitored during taping and playback on an Dage-MTI HR 
1000 monitor. It provides 1000 TV lines of resolution on a 30.5 cm (12 inch) flat face, 
cathode ray tube (CRT) display. The monitor offers selectable underscan, variable video 
enhance, switchable 110/120 VAC operation, and brightness and contrast controls. The 
monitor is operated at 110 VAC with normal underscan and the enhancement, bright¬ 
ness and contrast controls turned fully clockwise. The Panasonic VCR allows the vide¬ 
otape to be viewed one field (1/60 second) at a time. 
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VI. EXPERIMENTAL RESULTS 


This section presents surveys of measured total pressure, mean velocity, vortex cir¬ 
culation, mean vorticity and power spectra, as measured in a radial/spanwise plane lo¬ 
cated 120° from the start of curvature for Dean numbers from 35 to 450. Flow 
visualization results are also given for Dean numbers ranging from 50 to 425, and 
streamwise locations of 40°, 80°, and 120°. Data are collected at a total of 320 locations 
in the 6 = 120° radial/spanwise plane using a grid with eight radial (yjd) locations and 
40 spanwise (z/d) locations. The spanwise locations are referenced to the channel 
centerline, increasing from zjd = 4 to zjd = 8 and spanning 5.08 cm. The Kiel probe and 
five-hole pressure probe are used to make pressure measurements. Pressure measure¬ 
ments obtained with the five-hole probe are used to determine the components of ve¬ 
locity. The components of vorticity are determined from these velocity components. 

Power spectra results are collected at 6 = 112°, yjd = 0.5 and spanwise locations of 
9.84 cm or 10.16 cm from the channel side wail, which correspond to zjd = 12.25 and 
zjd = 12.0. 

Flow visualization results are collected at two spanwise locations (labeled inside and 
outside) at 6 = 40°. Inside refers to a portion of the radial/spanwise plane which extends 
from zjd = 3.94 to zjd = 7.87. Outside refers to a portion of the radial/spanwise plane 
which extends from z/d = 12.99 to zjd = 16.93. Each photograph in a particular se¬ 
quence shows the concave surface at the bottom and the convex surface at the top. The 
time interval between individual photographs in different sequences vary from 1/15 sec¬ 
ond to 1/60 second. 
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A. SURVEYS OF TIME-AVERAGED QUANTITIES 

Time-averaged results are presented in Figures 18 through 76. Of these, total pres¬ 
sure results from the Kiel probe are presented in Figures 18 through 35. Total pressure 
results from the five-hole probe are presented in Figures 36 through 40. Streamwise 
velocity contours are presented in Figures 41 through 44. Plots of velocity vectors in the 
curved channel flow are presented in Figures 45 through 47. Velocity perturbation 
contours are presented in Figures 48 through 51. Streamwise vorticity results are pre¬ 
sented in Figures 52 through 58. Streamwise circulation results are presented in Figures 
59 through 61. Radial vorticity results are presented in Figures 62 through 68. Spanwise 
vorticity results are lastly presented in Figures 69 through 76. 

1. TOTAL PRESSURE AND STREAMWISE VELOCITY 

Total pressure measurement results obtained using the Kiel probe are presented 
in Figures 18 through 35. These contour plots closely resemble results obtained with the 
five-hole pressure probe. Figure 35 shows linear increases of maximum and minimum 
values of total pressure with Dean number. 

Total pressure and streamwise velocity measurement results obtained using the 
miniature five-hole pressure probe are given in Figures 36 through 44. Figure 40 shows 
a nearly linear increase of five-hole total pressure measurements with increasing Dean 
number. Maximum and minimum refer to the maximum and minimum pressure differ¬ 
ences recorded for each Dean number. A comparison of the total pressure and 
streamwise velocity results indicates them to be qualitatively very similar for all Dean 
numbers investigated. Areas of low total pressure and low streamwise velocity develop 
near the concave wall as the Dean number increases from 35.03 to 61.92. Secondary 
flows generally move from the concave wall to the convex wall within these regions, 
which are present between pairs of counter-rotating vortices. They are referred to as 
upwash regions [Refs. 11,16], where upwash refers to the direction of flow with respect 
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to the concave surface. Plots of secondary flow velocity vectors which illustrate t^ese 
upwash regions are presented for Dean numbers of 86.69, 99.90, 149.8 in Figures 45 
through 47. At least three upwash regions are apparent in all three figures. As the Dean 
number increases above 149.8, Figures 36 through 43 show that the distance between the 
middle and right upwash regions decrease*. At De= 187.7, only two upwash regions are 
evident. Further ircreases in Dean number (Figures 39 and 44) result in a spanwise 
uniform flow. Flow visualization from the present study as well as from Ligrani and 
Niver [Ref. 11], show that vortex pairs exist at these higher Dean numbers, but are not 
apparent in time-averaged measurements because of considerate vortex motion in the 
radial and spanwise directions [Ref. 17]. 

2. VELOCITY PERTURBATION 

At small values of Dean number, the fluid motion in the curved channel is 
streamwise with no secondary flows. The velocity profile is nearly parabolic with a 
maximum which is shifted slightly closer to one wall (the concave surface)compared to 
flow in a straight channel. This is called curved channel Poiseuille flow (CCPF) [Refs. 
8,9]. To show mean velocity deviations from CCPF, values of CCPF are subtracted 
from the measured streamwise velocity. This difference is then normalized using the 
bulk mean velocity, U, for each Dean number studied to give the normalized velocity 
perturbation. Normalized velocity perturbation contours are presented in Figures 48 
through 51. Figure 48 shows that the flow changes from a nearly spanwise uniform one 
to one with spanwise periodic perturbations as the Dean number increases from 35 to 
about 100. The most s ; _nificant negative velocity perturbations occur at the same z/d 
locations as the upwash regions. Figures 50 and 51 show that the flow becomes 
spanwise uniform as the Dean number becomes greater than 250. 


23 




3. VORTICITY COMPONENTS 


The three components of vorticity are determined from the three components 
or' velocity using the computer program VORTICITY, which implements equations 
{10},(11}, and {12} from Chapter III. Figures 52 through 76 present the results of these 
calculations along with circulation magnitudes determined from the streamwise compo¬ 
nent of vorticity. The contour plots of this section show positive vorticity as continuous 
lines, and negative vorticity as dashed lines. In the radial;spanwise plane with the flow 
directed away from the observer, positive streamwise vorticity is counterclockwise and 
negative streamwise vorticity is clockwise. Graphs showing the variation of vorticity 
with Dean number for each of three vortex pairs apparent in the contour plots (left, 
middle, and right), appear at the end of the group of plots for that particular vorticity 
component. The minimum value of negative vorticity is presented as an absolute value 
on these plots. 

«. STREAMWISE VORTICITY 

Streamwise vorticity results are shown in Figures 52 through 58. 
Figure 52 shows that the flow does not display spanwise periodic areas of positive and 
negative vorticity across the span of the measurement plane for Dean numbers less than 
S6.69. At De = 99.9, the thrve regions containing both positive and negative vorticity 
are apparent. Each area of positive and negative vorticity evidences a pair of counter¬ 
rotating vortices, or a vortex pair As the Dean number increases, vorticity magnitudes 
in the pairs increase. The areas between the positive and negative vorticity regions in 
each pair are centered at the same zjd locations as the upwash regions. Figures 53 and 
54 show thr.t three pc.irs are replaced by two pairs as the Dean number increases from 
149.8 to 175.5. De= 175.5 results show a small region of negative vorticity located at 
yjd = 0.2, zjd ~ 6.3 which is very close to a i -gion of negative vorticity of a vortex pair 
located on the right-hand side of the measurement plane. These two areas of negative 
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vorticitv are replaced by one large region at De= 187.7 (Figure 54), which is centered 
at y/d = 0.4, z/d = 6.8. As these events occur, the spanwise locations of pairs also change 
somewhat with time. 

As the Dean number becomes greater than 250, areas of positive and neg¬ 
ative vorticitv are less apparent on the plots. The contour plots then show spanwise 
uniform time-averaged distributions of streamwise vorticity. Figures 56 through 58 show 
how the peak streamwise vorticity varies with Dean number. Plots are given for all three 
vortex pairs apparent in the contour plots. These are denoted as left, middle and right 
pairs according to their positions. Peak vorticity levels then generally increase with 
Dean number for all three cases. 

b. STREAMWISE CIRCULA TION 

Circulation values of streamwise vorticity, T, , are determined using the 
computer program CIRC1. This is done by summing the values of vorticity over a 
specified region, and then multiplying this result by the area corresponding to individual 
values of vorticity. Summing values of positive vorticity gives positive streamwise cir¬ 
culation and summing values of negative vorticity gives negative circulation. Figures 59 
through 61 present values of positive and negative streamwise circulation versus Dean 
number. Figure 59 values are given for the left vortex pair and show a general increase 
with Dean number up to 138.7. Above De= 138.7, the positive circulation decreases and 
then increases, while the negative circulation continues to increase as the Dean number 
increases further. Figure 60 shows an increase of both positive and negative circulation 
for the middle vortex pair with Dean number for values up to 127.4. At greater Dean 
numbers, values of positive and negative circulation tend to decrease. Values of positive 
and negative circulation for the right vortex pair in Figure 61 both generally increase 
with Dean number for the entire range of values shown. 
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c. RADIAL VORTICITY 


Radial vorticity distributions are presented in Figures 62 through 68. 
Spanwise periodicity is evident in results presented in Figures 62 and 63 for Dean num¬ 
bers ranging from 73.63 to 226.8. At Dean numbers greater than 250.3, the time- 
averaged flow is spanwise uniform. Figures 66 through 68 show the variation of peak 
radial vorticity with Dean number for the left, middle and right vortex pairs. For all 
three situations, the values of positive and negative vorticity generally increase until 
De= 149.8, and then tend to decrease as the Dean number increases further. 
d. SPANWISE VORTICITY 

Spanwise vorticity distributions are given in figures 69 through 76. Results 
in Figure 69 for Dean numbers of 35.03 and 51.36 show a flow field which is nearly 
spanwise uniform. As the Dean number increases, some portions of the contours of 
positive spanwise vorticity move toward the concave wall, while contours of negative 
spanwise vorticity move toward the convex wall. In Figure 73, the data of Figures 71 
and 72 are given in profile form. The time-averaged flow field becomes spanwise uni¬ 
form for values of Dean number greater than 250.3. Figures 74 through 76 show the 
variation of peak spanwise vorticity with Dean number. In all three figures, positive and 
negative vorticity values increase approximately linearly with Dean number for the left, 
middle, and right pairs for the ranges of values shown. 

B. SPECTRA OF HOT-WIRE SIGNAL MEASUREMENTS 

Spectral results are presented in Figures 81 through 94. The streamwise and radial 
locations of the probe for all measurements are 9 = 112° and yld = 0.5, respectively. The 
spanwise location for Figures 81 through 84 is 9.84 cm (3.875 in) from the channel wall 
or zjd = 12.25. The spanwise location for Figures 85 through 94 is 10.16 cm (4.0 in) 
from the channel wall or z/d = 12.0. 
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Figure 81 shows that the highest energy level for De = 40.3 is measured at 0.5-2.0 
Hz. At this Dean number, four other broad spectral peaks also exist at 5-6 Hz, 7-9 Hz, 
12-15 Hz and 17 19 Hz. As the Dean number increases up to 60.0, these peaks continue 
to be present as overall energy levels of the spectra increase. At De= 135.2, the highest 
energy level is measured at about 6 Hz. Only two other broad spectr 1 peaks exist at 
1-2 Hz and 12-14 Hz.. At De= 145.3, individual broad peaks are not so apparent. The 
time traces in Figure 82 show that the dominant frequency increases from about 1.5 
Hz (3 cycles /2 sec) for De = 40.3 to 5.5 Hz (11 cycles/2 sec) for De= 135.2. 

Figure 83 shows that the highest energy level is measured at 5-6 Hz for De= 140.6. 
Two peaks at 1.5-2.0 Hz and 11-12 Hz are similar to the ones seen in Figure 81. 
Figure 83 also show large broad spectral peaks at 60-100 Hz for Dean number from 
140.6 to 155.-1. This results because of twisting Dean vortex flow, described by Finlay 
et al. [Refs. 9,10]. Figure 84 presents the time traces associated with the previous figure. 
Here, peak-to-peak amplitudes increase with Dean number. 

Figure 85 show's a spectral survey for Dean numbers from 50.5 to 435.0. Individual 
spectra that make up this figure are also given in Figures 86 through 94. At De= 50.5, 
De= 100.3, and De= 150.1, a spectral peak at 60 Hz is present which is due to electronic 
noise. Figure 88 for De = 150.1 again shows a spectral peak at 80 Hz due to the pres¬ 
ence of twisting Dean vortex flow. 

C. FLOW VISUALIZATION 

Flow visualization results are presented in Figures 98 through 126. Table I gives 
Dean numbers and streamwise locations for the different figures. Flow visualization 
results are given for 6 = 40°, 80° , and 120°. Tw'o portions of the radial/spanwise (inside 
and outside ) are used. Inside refers to a radial/spanwise plane which extends from 
z/d = 3.94 to z/d = 7.87. Outside refers to a radial/spanw'ise plane which extends from 
zjd = 12.99 to z/d = 16.93. Time increases in each figure according to the frame numbers 
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to the left of each photograph. Each photograph shows a view with the concave surface 
on the bottom and the convex surface on the top. 


Table 1. DEAN NUMBER AND LOCATION OF FLOW VISUALIZATION SE¬ 
QUENCES 


DEAN 

NUMBER 

40° 

(inside) 

40° 

(outside) 

80° 

(inside) 

mmSMM 

El 

50 

Figure 98 

Figure 107 


Figure 119 

75 




Figure 120 

100 

Figure 99 



Figure 123 

125 




Figure 124 

150 

Figure 100 



Figure 125 

200 

Figure 101 

Figure 110 

Figure 116 

Figure 126 

250 

Figure 102 

Figure 111 

Figure 117 


300 

Figure 103 

9^901 

Figure 118 


350 

Figure 104 

Figure 113 



400 

Figure 105 

Figure 114 



450 

Figure 106 

Figure 115 




Photographic sequences of results obtained at 6 = 40°(inside) are now presented in 
Figures 98 through 106. Figure 98 (De = 50.9) shows evidence of flow in the spanwise 
direction. Dean vortex formation is not present as expected from the time-averaged 
contour plots. Figure 99 (De = 100.1) shows vortex pairs emerging from flow near the 
concave wall and disappearing into flow near the concave wall. Frames 7-17, 20-23, and 
31-33 show spanwise flow above the vortex pairs. Figure 100 (De= 150.7) shows addi¬ 
tional vortex pairs emerging from flow near the concave wall. Spanwise flow is clearly 
visible in frames 10-17. Figure 101 (De= 200.5) also shows a vortex pair emerging from 
flow near the concave wall in frames 3-5, and spanwise directed flow in frames 1-6. 
Frames 8-15 briefly show a spanwise periodic flow. Figure 102 (De= 250.8) shows the 
disappearance of a vortex pair due to its engulfment into the upwash of an adjacent 
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vortex pair in frames 22-32. Figure 103 (De= 300.3) shows a rocking motion of the 
vortex pairs. The poorly defined vortex pair in the middle of frame 1 moves slightly left, 
then right and then back to the left. It is more clearly evident in frame 30. Figure 104 
(De= 350.3) also shows this rocking motion. A vortex disappearance event is evident in 
frames 3-9. Here, two vortex pairs are present near the middle of frame 3, one slightly 
right and one slightly left. The positive vortex of the pair to the right then merges with 
the negative vortex of the pair to the left. Frame 9 then shows the outline of a single 
vortex pair which replaces the two pairs present previously. Figure 105 (De = 400.4) 
shows vortex movement which is too rapid for the video camera to capture without 
blurriness. Here the image is focused but appears slightly blurry. This blurriness also 
results because smoke particles do not follow all of the higher frequency vortex motions 
(i.e., twisting or turbulence). Figure 106 (De = 425.0) shows a higher degree of the 
chaotic vortex motion than seen at De = 400.4. A variety of appearance and disappear¬ 
ance events are evident along with other complex vortex motions. 

Photographic sequences of results obtained at d — 40 °{outside) are now presented in 
Figures 107 through 115. Figure 107 (De=50.9) shows approximately two cycles of 
spanwise flow where one cycle consists of flow from left to right and then back again. 
Figure 108 (De= 100.0) shows spanwise flow along with vortex pair formation. Frames 
22-34 show a vortex pair forming to the left of the spanwise midpoint of the photo¬ 
graphs. Figure 109 (De= 150.5) shows a more or less fully developed pairs as well as a 
vortex pair emerging from flow near the concave wall. The lifetime (over 2 seconds) of 
the more developed pair is unusually long for this Dean number. Typically, such pairs 
are present for less than one second. Figure 110 (De = 200.3) shows a rocking motion 
of the vortex pairs. Frames 7-15 show the emergence of a vortex pair from flow near the 
concave wall. Figure 111 (De = 250.6) also show's a rocking motion of the vortex pairs, 
as well as the disappearance (frames 3-14) of the vortex pair initially located in the 
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middle of frame i. Figure 112 (De= 300.7) shows vortex formation as well as spanwise 
flow. Frames 14-19 show spanwise flow around the exterior of a vortex pair. 
Figure 113 (De= 350.7) shows a vortex pair forming between two existing vortex pairs. 
The outline of this newly formed vortex pair is particularly evident in frames 31-34. 
Figure 114 (De = 400.9) shows vortex pair formation along with a violent rocking mo¬ 
tion in images which are somewhat blurred. Figure 115 (De = 425.3) shows vortex pair 
formation which is again present with fairly violent rocking motion. 

Photographic sequences of results obtained at 6 = 80 °{inside) are now presented in 
Figures 116 through 118. Figure 116 (De= 200.4) shows two interesting flow events. 
Frames 1-12 show one vortex pair splitting to form two vortex pairs. The outline of the 
initial single vortex pair is visible in frame 1. In frame 8, a downwash region is seen as 
a bright line at the spanwise centerline of the photograph, extending the full channel 
height. The downwash region is clearly visible in frame 12. The second flow event of 
interest is evident in frames 41-51, which show two adjacent vortex pairs joining to form 
one vortex pair. One vortex pair is located at the middle of frame 41, the other one is 
slightly to the left, and the upwash regions for each of the two vortex pairs are tilted 
toward each other. In frame 46 the two upwash regions are close enough to combine, 
resulting in one single vortex pair in frame 51. Figure 117 (De = 250.3) shows violent 
rocking of vortex pairs. In frames 6-11, the spanwise movement of the vortex pairs is 
so rapid that the images become blurry. Figure 118 (De= 300.3) shows similar behavior 
except that images are even more blurry. 

Photographic sequences of results obtained at 6= 120 °(inside) are now presented 
by referring to photographic sequences presented in Figures 119 through 126. 
Figure 119 (De=50.9) shows the formation of a vortex pair as it slowly drifts to the 
right. Frame 13 shows initial formation, just to the right of the spanwise centerline of 
the photograph. Figures 120 through 122 are photographic sequences of flows at 
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De= 75.5. Figure 120 shows two vortex pairs merging into one pair as the result of 
cancellation of individual vortices having opposite signs [Ref. 13]. Figure 121 shows two 
vortex pairs combining in frames 6-13 as a result of vortex pair engulfment [Ref. 13]. 
Frames 17-20 then show a new vortex pair emerging from flow near the concave wall. 
Frames 22-28 show that this newly formed pair then combines with the vortex pair di¬ 
rectly to the left, also because of pair engulfment. Figure 122 shows a vortex pair dis¬ 
appearing into flow near the concave wall [Ref. 13]. Figures 123 (De= 100.6) and 124 
(De= 125.5) show a gentle rocking of vortex pairs, while Figure 125 (De= 150.5) shows 
rocking vortex pairs which is higher in amplitude. Figure 126 (De= 200.8) show the 
emergence and disappearance of a vortex pair, which is located at the spanwise middle 
of the photograph in frame 8. 
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VII. SUMMARY AND CONCLUSIONS 


A transparent curved channel with a 40 to 1 aspect ratio w*as utilized to study the 
structural characteristics of Dean vortices at Dean numbers from 35 to 400. 

Mean velocity components and vorticity components are determined from time- 
averaged total pressure measurements measured using a miniature five-hole probe at a 
streamwise location of 6 — 120° . At low Dean numbers, (De= 35.03 to De= 61.92), 
these results show flow which is spanwise uniform. As the Dean number increases, evi¬ 
dence of vortex pairs is apparent in flow which appears spanwise periodic. The time- 
averaged flow becomes less spanwise periodic and more spanwise uniform as Dean 
numbers become greater than 250. 

Power spectra of voltage fluctuations from a hot-wire probe sensing velocity fluctu¬ 
ations were obtained at a streamwise location of 9 — 112°. Twisting Dean vortex flow 
is evidenced by spectral peaks measured for De= 140.6 to De= 155.4. The frequency 
of the principal spectral peak increases from about 70 Hz to 85 Hz as the Dean number 
increases over this range. 

Flow visualization results for streamwise locations of 6 — 40°, 80° , and 120°, show 
a variety of unsteady Dean vortex pair behavior, including, spanwise directed flow, 
vortex pair appearance, vortex pair disappearance, vortex pair merging, and vortex pair 
rocking. Spanwise flow is evident for De=50 and De= 100 for a streamwise location 
of 6 = 40°. Vortex appearance and disappearance events are also apparent, becoming 
more so at this location as the Dean number increases. Results at 8 = 80° for De= 200 
to De= 300 show vortex pairs w’ith rocking motion. At 6 = 120°, appearance and dis¬ 
appearance events are evident for Dean numbers of 50.9 to 200.8, and vortex pairs show 
rocking motion for Dean numbers of 100.6 to 150.5. As the streamwise location changes 
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from 0 = 40° to 6= 120°, appearance events, disappearance events, and vortex pair 
rocking occur at decreased Dean numbers. 


33 









SCREENS 


Z 



Figure 3. Schematic of Test Facility 
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Figure 4. Test Facility (Side View) 
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Figure 10. Pressure Probe Calibration Sled 
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Figure II. Wind Tunnel Used For Probe Calibration 









































Figure 12. Curved Channel Traversing Mechanism, Transducers, Probe and Mount 
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FIVE-HOLE PRESSURE PROBE 



S * E Z T 0 T— Z— E— t— 9- 

AdO 


Figure 13. C n vs. Yaw Angle, 3.9 m/s 
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Figure 18. Kiel Probe Total Pressure Contour, De= 51.81 
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Figure 19. Kiel Probe Total Pressure Contour, De = 77.47 
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Figure 20. Kiel Probe Total Pressure Contour, De= 103.6 
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Figure 22. Kiel Probe Total Pressure Contour, De= 149.6 
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Figure 25. Kiel Probe Total Pressure Contour, De = 226.8 
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Figure 27. Kiel Probe Total Pressure Contour, De = 274.5 
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Figure 28. Kiel Probe Total Pressure Contour, De = 302.0 
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Figure 29. Kiel Probe Total Pressure Contour, De = 327.5 
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Figure 30. Kiel Probe Total Pressure Contour, De= 350.7 
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Figure 37. Five-Hole Probe Total Pressure Contours, De= 111.8 to De= 187.7 
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VELOCITY VECTORS IN CURVED CHRNNEL FLOW 



in m m 

• CM 


P/A 

Figure 45. riot of Velocity Vectors for De = 86.69 
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Figure 46. Plot of Velocity Vectors for De = 99.9 
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figure 48. Velocity Perturbation Contours, De = 35.0 to De= Ili.8 
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Figure 49. Velocity Perturbation Contours, De= II 1.8 to De= 187.7 
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ignre 53. Streanmise Vorticity Contours, De= 111.8 to De= 187.7 
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Figure 54. Stream»ise Vorticity Contours, De= 187.7 to De = 275.4 
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Figure 55. Streamwise Vorticity Contours, De= 275.4 to De-428.6 
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Figure 56. Peak Streannvise Vorticity vs. Dean Number, Left Vortex Pair 
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Figure 62. Radial Vorticity Contours, De=35.0 to De= 111.8 
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Figure 63. Radial Vorticity Contours, De= 111.8 to De= 187.7 
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Figure 64. Radial Vorticity Contours, De= 187.7 to De = 275.4 
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Figure 65. Radial Vorticity Contours, De = 275.4 to De = 428.6 
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Figure 66. Peak Radial Vorticity vs. Dean Number, Left Vortex Pair 
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Figure 67. Peak Radial Vorticity vs. Dean Number, Middle Vortex Pair 
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Figure 68. Peak Radial Vorticity vs. Dean Number, Right Vortex Pair 
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Figure 73. Spanwise Vorticity vs. Radial Location 
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Figure 75. Peak Spamvise Vorticity vs. Dean Number, Middle Vortex Pair 
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Tigure 76. Peak Spanwise Yorticity vs. Dean Number, Right Vortex Pair 
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Figure 78. Overview of the HP 3562A Signal Analyzer Measurement Process 
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Figure 79. Typical Settings Used on the HP 3562A Signal Analyzer 
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Figure 83. Conventional Hot-wire Probe Spectra, De= 140.6 to De= 155.4 
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Figure 84. Conventional Hot-» 
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Figure 85. Conventional Hot-wire Probe Spectra, De= 50.5 to De = 435.0 
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Figure 86. Conventional Hot-wire Probe Spectrum, De = 50.5 
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Figure 90. Conventional Hot-wire Probe Spectrum, De = 250.0 
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Figure 92. Conventional Ilot-uire Probe Spectrum, De= 350.3 
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Figure 95. Flow Visualization Radial/Spanwise Plane 
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Cooled Smoke to Channel Inlet 




96. Schematic of Smoke Generator and Smoke Rake 
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Figure 104. Flow Visualization, De = 350.3, 6 = 40°(inside), 1/60 Second Intervals 
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Figure 120. 



Flow Visualization, De = 75.5, 0 = 120°(inside), 1/30 Second Intervals 
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Figure 122. 



Flow Visualization, De = 75.5, 0 = I20°(inside), 1/30 Second Intervals 


154 

























155 















































Figure 125. Flow Visualization, De= 150.5, 0 = 120°(msi<fe), 1/30 Second Intervals 
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APPENDIX B. UNCERTAINTY ANALYSIS 


The uncertainty analysis presented below is taken from Baun [Ref. 17 ]. The present 
study deals with flow situations similar to those explored by Baun, and utilizes the same 
facilities and equipment. 

A. DEAN NUMBER UNCERTAINTY 

The uncertainty associated with the calculation of the Dean number is derived from 
the uncertainty in determining the pressure drop across the ASME orifice plate. Values 
of uncertainty for De= 150 are listed below. 

1. CELESCO TRANSDUCER CALIBRATION UNCERTAINTY 
3P„„, = ± 0.002 in. H 2 0 

SE = ± r .0007 Volts 
5C = ±0.0012 in. H 2 OfVolt 

2. PRESSURE DROP UNCERTAINTY 
3E = + 0.047 Volts 

<5AP or = ± 0.005 in. H } 0 

3. MASS FLOW RATE UNCERTAINTY 
SA or = ± 0.012 in 1 

6K=± 0.001 

3Y= ±0.02 

dp — ± 0.002 Ibmlft 2 

dm = ± 0.0004 Ibml sec (95% confidence level) 

3 A c/I = ± 0.0052 in 2 

3De = ± 4.4 (95% confidence level) 






B. FIVE-HOLE PROBE VELOCITY MEASUREMENT UNCERTAINTY 

The uncertainties associated with the Five-hole probe velocity components are de¬ 
rived from the uncertainty in probe calibration as well as the uncertainty with the pres¬ 
sures measured with the five-hole probe. Values of uncertainty for De= 100.3, 
ytd— 0.5, z/d— 6.0 are listed below. The probe calibration uncertainty is for 
yaw = —10°, pitch = 0°. 

1. VALIDYNE TRANSDUCER CALIBRATION UNCERTAINTY 

= ± 0.001 HI. H 2 0 
SE = ± 0.0007 Volts 

6C = ± 0.000072 in. H 2 OIVolt (bias error is disregarded since it will cancel out 
when pressure differentials are calculated) 

2. PROBE CALIBRATION UNCERTAINTY 
SE = ± 0.0004 Volts 

<5P, = ± 0.0006 in. H 2 0 
dP 2 = ± 0.0001 in. H 2 0 
SP 2 = ± 0.0005 in. H 2 0 
6P t = ± 0.0002 in. H 2 0 
6P S = ± 0.0003 in. H 2 0 
5P S = ± 0.00004 in. H 2 0 
6P,= ± 0.0007 in. H 2 0 
5P = ± 0.00015 in. H 2 0 
<5C„ = ± 0.025 
SC„ = ± 0.011 
6C,„ = ± 0.045 
dK y = ± 0.004 degrees "■ 
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SK P = ± 0.004 degrees _l 
SK„ = + 0.003 degrees' 1 

3. CURVED CHANNEL MEASUREMENT UNCERTAINTY 

SE = ± 0.0007 Volts 
SC = ± 0.000072 in. H 2 0/Volt 
SP i = ± 0.0002 in. Ii 2 0 
SP 2 = ± 0.0002 in. H 2 0 
SP 3 = ± 0.0002 in. H 2 0 
SPi = ± 0.0002 in. H 2 0 
SP s = ± 0.0002 in. H 2 0 
SP = ±0.0001 in. H 2 0 
SC„ = ± 0.134 
SC ep = ± 0.141 
S Yaw = ± 1.33° 

S Pitch = ± 1.4° 

S(P,-P,) = ± 0.0004 in. H 2 0 
<5U=± 0.063 m/s 
S U, = ± 0.062 mjs 
SU, — ± 0.062 m/s 
SU, = ± 0.030 m/s 
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APPENDIX C. SOFTWARE DIRECTORY 


This appendix lists the various programs used in this study. Accompanying each 
program listed, is a general summary of what the program does, user inputs, program 
outputs and additional features. Programs are listed alphabetically by usage: Kiel probe, 
five-hole pressure probe, and plotting routines. All programs are written/modified in 
BASIC 4.0 for use on the HP 300 Model 9000 computer.!Ref. 17] 

A. KIEL PROBE PROGRAMS 

KIEL15: This program is used to control Kiel probe positioning and data acquisition 
in the curved channel. It interfaces with the MITAS controller to control the probe 
position and is similar to the program FIVHOLE but obtains data for a single pressure. 
The i5 in the program name refers to the i.5 inch orifice plate used between the outlet 
and blov/er plenums. 
user input 

• estimated Dean number. 

• number of spanwise and radial locations at which data is to be taken. 

• spanwise and radial resolution. 

• initial position of the probe. 

• number of samples to be taken per channel. 

• number of standard deviations. 

• number of consecutive runs. 

• transducer calibration coefficient. 

• ambient pressure. 

• time to delay the start of the experiment. 


program output 

« all information input by the user. 

• data file of raw pressures: y, z, P ]t De. 

• average Dean number. 
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additional features 

• reads temperature and calculates ambient density and Dean number for each 
measurement location. 

• allows user to calibrate transducer or input the calibration coefficient. 

B. FIVE-HOLE PRESSURE PROBE PROGRAMS 

ALIGN: This program aligns the five-hole pressure probe in the curved channel at zero 
yaw and zero pitch. It also interfaces with the MITAS controller to position the probe 
at a desired location and repositions the probe to the original position upon completion 
of the program. 
user input 

• number of samples per channel. 

• number of standard deviations. 

• ambient pressure. 

• transducer calibration coefficients. 

• original position of the probe. 

• desired position of the probe. 

program output 

• mean velocitv. 

• A. Pi, P>. P*\Py 

• C„, C„. 

additional features 

• allows the user to realign the probe between measurement sequences. 

• allows user to calibrate transducers or input their calibration coefficients. 

CCPF1: This program calculates the Curved Channel Poiseuille Flow (CCPF) at a given 
Dean number. 
user input 

• data file of velocities as output by the program VFLOCITY1. 
program output 

• data file of streamvrise velocities: y, z, UCCPF, U,— UCCPF . 

additional features 

• none. 
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C1RC1: This program calculates the peak vorticity levels and circulation values from a 
particular set of vorticity data, using user defined boundaries. 
user input 

• data file of vorticity as output by the program VORTICITY. 

• threshold limits, as a percentage, for omitting data points when circulation is 
calculated. For example, a 20 percent threshold will cause the program to omit 
all vorticity values less than 20 percent of peak vorticity. 

• spanwise (z/d) locations that define boundaries of the areas of interest. 


program output 

• hard copy and or binary data file of y/d and z/d limits, peak vorticity and circu¬ 
lation for each vortex pair defined by the user. 

additional features 

• none. 


DEAN 15: This program samples the orifice pressure drop and calculates the Dean 
number as discussed in Chapter II. The 15 in the program name refers to the 1.5 inch 
orifice plate used between the outlet and blower plenums. 
user input 

• ambient pressure. 

• transducer calibration coefficient. 


program output 

• transducer calibration coefficient if calculated. 

• expansion coefficient, Y. 

• How coefficient, K. 

• orifice pressure drop. 

• mass flow rate. 

• channel Reynolds number. 

• Dean number. 


additional features 

• allows user to adjust the system throttle valve (flow rate) and then calculate an¬ 
other Dean number. 

• allows user to calibrate transducer or input the calibration coefficient. 
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FIVHOLE: This program is used to control the five-hole probe positioning and data 
acquisition in the curved channel. It interfaces with the M1TAS controller to control 
the probe position. 
user input 

• estimated Dean number. 

• orifice size. 

• number of spanwise and radial locations at which data is to be taken. 

• spanwise and radial resolution. 

• initial position of the probe. 

• time required for pressure reading to stabilize. 

• number of samples to be taken per channel. 

• number of standard deviations. 

• transducer calibration coefficient's. 

• ambient pressure. 

• time to delay the start of the experiment. 


program output 

• all information input by the user. _ 

• data file of raw pressures: y, z, P, , P 2 , P 3 , P 4 , P s , P , De. 

• data file of uncorrected pressure coefficients: y, z, C py , C pp . 

• average Dean number. 


additional features 

• reads temperature and calculates ambient density and Dean number for each 
measurement location. 

• allows user to calibrate transducers or input their calibration coefficients. 

ORIENT: This program is used to orient the five-hole probe in the calibration wind 
tunnel before calibration begins. 
user input 

• number of samples to be taken per channel. 

• number of standard deviations. 

• transducer calibration coefficients. 


program output 

* mean velocity. 

* P|f P(> Pj> P«» P5 • 

* C„, c„. 
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additional features 

• allows user to continue through the program adjusting the probe calibration sled 
until satisfactory C py and C„ are attained. 

• allows user to calibrate transducers or input their calibration coefficients. 

PADJUST: This program implements the spatial resolution correction as described by 
Ligrani et al. [Ref. 16]. Chapter III of this study briefly describes this correction proce¬ 
dure. 

user input 

• distance between port 1 and the center of the other pressure ports, ylen and zlen. 
For the miniature five-hole pressure probe used in this study, both ylen and zlen 
are equal to 0.018 in. 

• data file of the raw pressure data as output by the program FIVHOLE. 


program output 

• file of corrected pressure coefficients: y, z, C py , C pp , P„ P, De. 


additional features 

• can be used with a survey of one spanwise location, i.e., a profile. 

PROCAL: This program is used for the calibration of the five-hole pressure probe in 
the calibration wind tunnel. 
user input 

• number of yaw points. 

• number of pitch points. 

• number of samples to be taken per channel, 

• number of standard deviations. 

• transducer calibration coefficients. 

• yaw and pitch angle resolution. 


program output 

• data file of calibration coefficients: pitch, yaw, C„, C pp , C p „ C p „ C p „. 

• data file of pressures: pitch, yaw, P u P 3 , P,, P 4 , P } . 

• P P 

1 1 r 


additional features 

• allows user to calibrate transducers or input their calibration coefficients. 

• allows user to vary yaw at constant pitch or vary pitch at constant yaw. 
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• can accommodate up to seven transducers. 

• allows for simultaneous printout of data if desired. 

VELOCITY1: This program computes the three mean velocity components from pres¬ 
sure data that has been corrected for spatial resolution. 
user input 

• data files of calibration coefficients and pressure data as output from program 
PROCAL. 

• number of yaw points. 

• number of pitch points. 

'■ • data file of pressure coefficients as output by program PADJUST. 

• values of y and z direction downwash velocity correction, dely and delz. In this 
study, a value of 0.011 is used for both dely and delz. 

program output 

• data file of pressures and velocities for each measurement location: y,z, P„ 

U, u, , U r , L\. 

additional features 

• can be used with a survey of one spanwise location, i.e., a profile. 

VORTICITY: This program calculates the three mean components of vorticity at each 
measurement location. Methods used for calculating vorticity are presented in Chapter 
III. 

user input 

• data file of velocities as output by the program VELOCITY1. 
program output 

, • data file of vorticity as follows: y, z , a> 9 , u)„ a >,. 

* additional features 

• can be used with a survey of one spanwise location, i.e., a profile. 

C. DATA PROCESSING 

CONTOUR: This program is a general contour plotting routine that has been adapted 
for specific plotting purposes. It plots contour lines between predetermined levels. 
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user input 


• number of data files to be averaged. 

• piot option: average values, average of deviations from mean value, or average 
of standard deviations from mean value. 

• number of contour levels to plot. 


program output 

• contour plot and table of levels. 

• maximum and minimum values. 

• average Dean number. 


additional features 

• none. 

PROFILE1: This program is a modified version of PLTRAV used by Baun [Ref. 17]. 
It has been modified to print out specific parameters for five-hole probe surveys con¬ 
ducted at one spamvise location, i.e., a profile. 
user input 

• data file of pressures and velocities as output by the program VELOCITY1. 

• column number of dependent variable. 

program output 

• spamvise location. 

• specified parameters with corresponding radial position. 

additional features 

• none. 

3DPLOT1: This program plots up to seven contour plots from curved channel meas¬ 
urements in a three dimensional format. An example of this format is presented in 
Figure 36. 
user input 

• tvpe of data plot desired: vorticity, streamwise velocity, total pressure, or 
U-UCCPF. 

• appropriate data files. 
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program output 

• three-dimensional presentation of specified contour plots. 

• table of levels. 

additional features 

• none. 


VECTOR!: This program plots secondary flow vectors present in the flow field. 
user input 

• data file of velocities as output by the program VELOCITYl. 


program output 

• plot of secondary flow vectors. 

• average Dean number. 

• actual length of largest vector and its corresponding velocity. 


additional features 

• allows user to zero velocity vectors present in the bottom row. 

• allows user to zero velocity vectors present at specific locations. 
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